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1. Introduction 

Compact stars are dense enough that they could contain novel phases of dense, strongly- 
interacting matter in their interior. In particular quark matter presents a challenging possibility 
[1]. To probe the physics of the core of a compact star is complicated since all observed electro- 
magnetic radiation originates exclusively from its surface or its magnetosphere and the physics of 
the crust in between these regions is complicated and poorly constrained. The only direct way to 
"see inside" the star is to use the same methods by which we know about the interior of the earth 
or sun, namely to study seismic oscillation modes. In contrast to the latter cases, these oscillations 
cannot be observed directly in a compact stai\ Yet, in case they have a non-vanishing quadrupole 
moment, like for instance r-modes, they couple to gravitational waves that carry away angular mo- 
mentum and can quickly spin down the star. For pulsars the rotation is directly observed via the 
pulsed emission which is among the most precise data in physics. Oscillations are damped by dissi- 
pative mechanisms, which relates the observable pulsar timing data on stars with large frequencies 
to the microscopic properties of the matter in the star's interior. 

The observed spin frequencies and spindown rates of pulsars show a striking age dependence. 
Young stars in particular feature spin frequencies that are more than an order of magnitude below 
their theoretical Kepler limit. Unstable r-mode oscillations [2, 3] have been suggested as an expla- 
nation for this surprising finding. Using semi-analytic expressions that allow to study uncertainties 
in the microscopic as well as macroscopic description, we have recently shown that r-modes of 
neutron stai^s can indeed provide a quantitative explanation for the low frequencies of young pul- 
sars [4]. In contrast to the final frequency the spindown time depends strongly on the amplitude 
at which the growth of the unstable oscillation is saturated by non-linear, amplitude-dependent 
damping processes and to obtain a quick spindown requires lai^ge saturation amplitudes. 

Old x-ray and radio pulsars, in contrast, are or have been spun up by accretion in binary sys- 
tems and can reach large rotation frequencies. They are very stable and their observed spindown 
rates are orders of magnitude smaller than those of young stars. With a saturation amplitude re- 
quired to explain the low frequencies of young stars, millisecond pulsars would spin out of the 
instability region on time scales that are significantly smaller than their age. Here we extend our 
recent analysis [4] and discuss different scenarios to make the r-mode mechanism for the spindown 
of young pulsars consistent with the observed data on old pulsars. One possibility is a temperature 
dependence of the r-mode saturation amplitude, but we show that even a strong temperature depen- 
dence does not allow to make the different observations compatible. We propose another scenario 
based on deconfined quai^k matter that has the potential to explain the timing data on both young 
and old pulsars. 

2. R-mode spindown of neutron stars 

In this section we discuss semi-analytic solutions for the spindown evolution of pulsars that are 
discussed in more detail in [4]. The basis of this approach is the observation that the microphysical 
material properties, like the shear and bulk viscosity rj and ^, the specific heat cy and the neutrino 
emissivity e that describe a given form of dense matter, depend on the macroscopic parameters 
relevant for the r-mode evolution, namely the star's temperature T and rotation frequency Q., over 
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the relevant parameter ranges via power laws 

71=7)7" , C~ — ' Cy=Cyr" , £ « g^^ (2.1) 

where the subthermal, low temperature Umit for the bulk viscosity has been employed which holds 
for fast spinning neutron stars at 7 ^ lO'" K [5] and the relevant oscillation frequency in the 
rotating frame is given for an r-mode by w 2/3 In observations pulsars appear as point sources 
so that only quantities that are averaged over the star are observable. The simple analytic parameter 
dependence allows to write the macroscopic quantities that are relevant for the r-mode evolution - 
given by the power radiated in gravitational waves Pg , the dissipated power Ps and Pg due to shear 
and bulk viscosity, the specific heat of the star Cy , the total neutrino luminosity Ly and the moment 
of inertia 7 - in a way that the complete microphysics as well as the considered star configuration 
is encoded in coiTcsponding dimensionless constants /, S, V, Cy, L and / (see [4] for more details). 
Qualitative different forms of matter can have very different power law exponents. For a standard 
neutron star the sheai^ viscosity is dominated by leptonic scattering processes [6] with (negative) 
power law exponent a = 5/3 and the neutrino emissivity is due to modified Urea reactions [7] with 
exponent = 8. 

The star evolution for the temperature T, the rotation frequency D. as well as the dimensionless 
global r-mode amplitude a are obtained from general energy and angular momentum conservation 
laws [8, 9]. An r-mode is unstable to the emission of gravitational waves when the growth is faster 
than the viscous damping 1/% > l/ty + I/T5, where T,- = Pi/ {2Em) are the corresponding charac- 
teristic time scales given in terms of the mode energy £",„. Then the amplitude grows exponentially 
on time scales of the order of seconds and this growth is eventually stopped by a non-linear satu- 
ration mechanism. There are several proposed mechanisms [10, 11, 12, 13], but at this point it is 
not settled which mechanism dominates. We will merely assume that such a mechanism exists that 
saturates the mode at a finite amplitude ttsat- In general the saturation amphtude is a function of 
temperature and frequency which we will model here by a simple power law dependence 

a^^^{T,Q) = a,^{r^Qj . (2.2) 

At saturation and for physical amplitudes otsat < 1 the residual evolution equations simplify to 

dO. IQaLO. dT 1 , 

Note that the required dissipation to saturate the r-mode, which is entirely determined by grav- 
itational physics, eventually ends up as heat and appears on the right hand side of the thermal 
evolution equation. 

In [4] it is shown that in general the thermal evolution of a compact star within the r-mode 
instability region is always faster than the r-mode spindown. The spindown of a compact stai^ 
follows thereby a steady-state curve where heating and cooling exactly balance. At temperatures 
that are high enough that neutrino cooling from the core dominates photon emission from the 
surface this curve is given for the dominant m = 2 r-mode and a general form of dense matter by 
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A similar semi-analytic expression has been previously given for the low temperature bound- 
ary of the instability region [14], where shear viscosity dominates the damping. The intersection 
of the two curves gives then the endpoint of the r-mode spindown evolution, where the star leaves 
the unstable region and the r-mode is damped. For a standard neutron star {NS) and assuming a 
constant saturation amplitude (/3 = 7=0) we find 



where the parameters with the suffix "fid" are the corresponding values for a fiducial 1 AMq neutron 
star with an APR equation of state [15] (see [14] for more details on the neutron star model). Note 
the extreme insensitivity of these expressions to the parameters characterizing the poorly known 
microphysical details like the shear viscosity {§) and the neutrino emissivity (L), as well as to the 
saturation amplitude ttsat- Fig- 1 compares the semi-analytic expressions to a numeric solution for 
the evolution of a young 1.4M0 neutron star for different constant saturation amplitudes (j8 = 7 = 
0). As can be seen the agreement is very good. The computation of the spindown time requires a 
solution of the spindown equation eq. (2.3), which gives in case of a neutron stai^ 

(2.7) 

In contrast to the endpoint of the evolution the spindown time depends strongly on the saturation 
amplitude. Whereas the spindown is fast for a^ai ^ 10^^ it becomes longer than the age of observed 
young pulsars for ttsat ^ 10^^. 

To compare a standard, static r-mode instability analysis to pulsars requires a temperature 
measurement which is only available for x-ray pulsars that currently accrete. Such an estimate 
involves many uncertainties since it requires to connect the core temperature to the observed surface 
temperature and thereby to model the heat transport in the complicated neutron star crust [16]. 
However, using the semi-analytic results for the spindown evolution it is possible to connect the 
results of the dynamic r-mode evolution to precise pulsar timing data. This requires the spindown 
rate eq. (2.3) which is given along the spindown curve eq. (2.4) by 



dt V3^52 

The comparison of the r-mode spindown evolution to pulsar timing data taken from the ATNF 
catalog [17] is shown in fig. 2. Here we extend the coiTcsponding plot given in [4] to include old 
pulsars which feature much lower spindown rates. The different lines show the r-mode spindown 
evolution of a 1.4M0 neutron star for different saturation amplitudes ranging from ttgat = 10^^ 
(bottom) to ttsat = 1 (top). Their endpoints map out the boundary of the instability region denoted 
by the corresponding line with an error band estimated from the uncertainties of the different pa- 
rameters in eq. (2.6) (see [4] for more details). Let us first consider young stars that have large 
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Figure 1: The spindown evolution of a young lAMp, neutron star with an APR equation of state [15] 
in temperature-angular velocity space. The boundary of the instability region of the fundamental m — 2 
r-mode is shown by the dotted curve. The dashed curves (which are mostly hidden underneath the solid 
curves) represent the steady state eq. (2.4) where heating equals coohng and are given for different constant 
(j3 =7 = 0) r-mode amplitudes ranging from ttsat = 10^^ (left) to Usat = 1 (right). The dot-dashed curve 
shows the steady state curve for j3 = 3.5 and asat = 10^^ at the boundary (see text). The solid lines show 
the numerical solution of the evolution equations for two fiducial initial spin frequencies Q. = 0.8 H^- and 
Q. = 0.2 Q.K, where Q.k is the Kepler frequency. As can be seen, after an initial cooling phase the evolution 
simply merges on the appropriate steady state curve and follows it to the edge of the instability region. The 
dots denote the semi-analytic results for the endpoints of the spindown evolution where the star leaves the 
instability region, eqs. (2.5) and (2.6). 

spindown rates and appear in the upper part of the plot. As can be seen all young stars are likely 
already outside of the instability region of a neutron star with standard damping mechanisms. Cor- 
respondingly, the r-mode mechanism in neutron stars can provide a quantitative explanation of the 
low spin frequencies of young pulsars. This could be a mere coincidence, but the quantitative 
agreement with the precise r-mode prediction is quite intriguing. In contrast, stars involving other 
forms of dense matter with enhanced viscous damping or other dissipative mechanisms in the crust 
have smaller r-mode instability regions (in a standard T-^l diagram) [18] and thereby generally 
cannot spin down pulsars to similarly low frequencies. 

As seen in fig. 2, most old stars spin at frequencies that are too low for r-modes to be relevant 
since other mechanisms, like magnetic braking, will spin these stars down to even lower frequencies 
(stars clustering around a diagonal line in fig. 2). Yet, when a star is in a binary system it can be 
spun up again by accretion from a low mass companion in a quasi-stable orbit over long time scales 
and reach nearly kHz frequencies (stars clustering along a horizontal branch in the lower part of 
fig. 2). These millisecond pulsars are clearly within the instability region of a standard neutron star 
and their spindown rates are many orders of magnitude lower than those of young stars discussed 
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Figure 2: The spindown evolution of a neutron star compared to observed pulsar data from the ATNF catalog 
[17]. Circles represent stars that spin down whereas squares represent stars that spin up. The solid lines show 
the evolution of the lAM^ star for different constant saturation amplitudes ()3 =7=0) ranging from Usat ~ 1 
(top) to ttsat = 10^^ (bottom). The dot-dashed line represents the steady-state for j3 =3.5 and ttsat = 10^^ 
at the boundary. The vertical dashed line to the left shows the lower frequency boundary of the instability 
region. The steep solid line is formed by the endpoints of the evolution for different values of a and thereby 
represents the boundary of the instability region in /-/-space, whereas the gray band reflects the error due 
to the uncertainty in the underlying parameters. The dotted line segments show the current evolution for two 
stars for which a reliable braking index is available [20, 21]. 



so far. Yet, large amplitude r-modes would have spun these stars out of the instability region on 
time scales much smaller than their current age of up to billions of years. To circumvent this, 
r-modes in neutron stars would have to saturate at very low (constant) amplitudes ttoid < 10^' in 
order to be compatible with the pulsars data. Similar bounds have e.g. recently been obtained from 
x-ray pulsars [19]. Since pulsars also spin down due to other mechanisms like magnetic dipole 
emission, those other mechanisms might then explain the scattered data e.g. by different magnetic 
fields. However, none of the currently proposed saturation mechanisms [10, 1 1, 12, 13] can saturate 
r-modes at such low saturation amplitudes. 
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3. Compatibility of the spin data of young and old stars 

Having seen that the r-mode mechanism in neutron stars can provide a quantitative explanation 
for the spin frequencies of young pulsars for sufficiently large saturation amplitudes, but would 
have problems to explain the comprehensive data on old pulsars since these would be well within 
the instability region, it is interesting to study if there is a consistent scenario that can describe both 
the r-mode spindown of young pulsars and the observed data on old millisecond pulsars. 

According to the results for constant saturation amplitudes in fig. 2, a possibility to explain 
this finding might be a strong temperature dependence of the saturation amplitude ttsat so that the 
spindown rate is large for hot young stars but many orders of magnitude smaller for cold old stars to 
meet the strict observational bounds. A qualitatively similar temperature dependence seems to be 
found in mode coupling models [10, 11] which present probably the most promising candidate for 
the r-mode saturation. Within such an explanation of the data the r-mode should be unimportant for 
the spindown evolution of old stars and other spindown mechanisms like magnetic braking would 
dominate. In order to model such a temperature dependence of the saturation amplitude we need 
large amplitudes 0.01 < ttyoung ^0.1 at core temperatures 10^ K < Tyoung ^ lO'^K relevant for 
young stars and very small values 10^^ < ttoid ^ 10^' at core temperatures 5 • 10^ K < Toid < 5 • 
10^ K realized in old stars. I.e. a change ttyoung/ cCoid ~ 10^ over a temperature range ryoung/^oid ~ 
50 which requires a large exponent jS 3.5 over this range. As can be seen from eq. (2.4) the 
steady-state curve along which the star spins down is thereby drastically changed from ~ T 
for a constant saturation amplitude (in case of neutron star) to the much weaker dependence D, ~ 
r^/^. Clearly such a strong temperature dependence of ttgat would need to be restricted to the 
above temperature interval and could not extend up to the boundary of the instability region in 
fig. 1 since this would give unphysically large amplitudes ttsat ^ 1 at temperatures T > lO^'^K. 
Assuming that the r-mode amplitude is constant above the considered temperature range, the star 
would quickly spin down along one of the large amplitude curves in fig. 1 and then merge into 
a more flat segment as shown by the dot-dashed curve for j8 = 3.5 along which it would exit the 
instability region after long times. Now let us estimate how the corresponding path would look 
in fig. 2. According to eq. (2.8) the spindown rate would have an extremely strong dependence 
on the frequency d£l/dt ~ Q.^^ for such a temperature dependent saturation amplitude leading to 
a nearly vertical segment shown by the dot-dashed curve in fig. 2. Clearly such a spindown path 
would be completely inconsistent with the low spindown rates of millisecond pulsai^s. Note that we 
assumed here neutrino cooling, but in case photon emission from the surface dominates for some of 
these sources this would not qualitatively alter the demonstrated effect. In other words r-modes in 
a neutron star with standard damping mechanisms cannot simultaneously explain both the low spin 
frequencies of young pulsars and the high spin frequencies of many old x-ray and radio pulsars. 

There is a more exotic but also more interesting possibility. During its initial spindown a young 
neutron star contracts since centrifugal forces become weaker. This increases the density and could 
at the con^esponding critical value transform hadronic matter in its interior into exotic forms of 
matter like hyperons or some form of quark matter, whereby a hybrid stai^ is formed. Such forms of 
matter have a very different bulk viscosity due to non-leptonic strangeness changing reactions with 
a resonant maximum at significantly lower temperature [5] , imposing a stability window where the 
star is stable against gravitational wave emission by r-modes up to large frequencies [14]. However, 
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in general the subsequent spin-up phase of recycled pulsars in binaries would reverse the process 
(transforming the star back to neutron matter) and could therefore not simultaneously explain both 
the spindown of young and old stars. Yet, when the strange matter hypothesis holds and strange 
quark matter is the true ground state of matter [22], once a quark core is formed it would be 
absolutely stable and could not be transformed back. In fact it would even grow and transform 
most of the star into strange matter. As has recently been shown [23], including important long- 
range strong interaction corrections to the bulk viscosity, the observed x-ray pulsars are outside of 
the instability region of a compact star that contains quark matter. Correspondingly, such a quark 
star could explain the large spin frequencies of observed x-ray pulsars, and as will be discussed in 
a forthcoming work it also has the potential to explain the data on radio pulsars. Neutron stars that 
are transformed into quark stars at a later stage of their evolution could thereby provide a universal 
explanation of the observed pulsar timing data. 
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